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As described in previous reports, this equation is constrained to the liquid-vapor coexistence boundary. 
Integration of the "thermodynamic equation of state" along isotherms involves a term for the compressibility factor 
of saturated vapor in the limit of low densities, which in previous work was not adequately defined. The present, 
brief report assumes familiarity with previous work, and presents a solution for the problem, which is to describe 
saturated-vapor densities in terms of the compressibility factor for saturated vapor, utilizing a given, vapor- 
pressure equation. 
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Symbols and Units 

Subscripts c and t refer to critical and liquid triple points. 

Subscript cr refers to liquid-vapor coexistence. 

d density, mol/L 

E the internal energy, J/mol 

J the joule, 1 N-m 

L the liter, 10~ 3 m 3 

mol 28.054 grams of ethylene (C 12 scale) 

P pressure in bars, 1 bar = 10 5 /V/ra 2 , (1 atm = 

1.01325 bar) 
P a (T) vapor pressure, bar 

P*(P) PJUp)] for equation (1) 

R the gas constant, 8.31434 {J/mo\)/K, 0.0831434 

(bar-L/mol)/A: 
R* (0. 083 1434) -d c , bar/K 

p d/d c , reduced density 

T temperature, K 

T a (p) liquid-vapor coexistence temperature 

1 . Introduction 



Adaptations of this type of equation of state have been 
described in numerous reports [1-8], J which must serve the 
interested reader as background, to avoid repetition here. 
The form of this isochoric equation is 

P - P a {p) =pR* ■ [T - T tT (p) ] + P 2 R*T C ■ F(p,T), (1) 

in which, for any given density (isochore), the coexistence 



temperature, T (T (p), is obtained by iteration from equations 
for the orthobaric densities. This is used to obtain the vapor 
pressure, P lT \T„(p) ], as a function of density. The function 
F(p,T) contains terms roughly analogous to T and 1/T 2 of the 
well-known, Beattie-Bridgeman equation of state, each with 
a density-dependent coefficient. Symbols and units for this 
report are given in a list. Fixed-points, used for some 
examples on ethylene, are given in table 1. 

Table 1 . Fixed-Point Values Used for Ethylene 





Triple Point 


Critical Point 


Temperature, K 


103.97 


282.35 


Pressure, bar 


1.2095- 10" 3 


50.4234 


Density mol/L 






Vapor 


1.3994- 1(T 4 


7.60 


Liquid 


23.3590 


7.60 



2. The Internal Energy 

Placing eq (1) in the thermodynamic equation of state 
yields 



A£ 



= f[(Z a 

Jo 



(p)- l)'RT a (p)/p + ...]-djp, (2) 



* This work was supported I by the American Gas Association, Inc., 1515 Wilson pressures p (f) and for the Saturated Vapor densities 
Boulevard, Arlington, VA 22209. 

1 Figures in brackets refer to literature references at the end of this paper. 



where E is the internal energy. 

Our problem is in the definition and value of 

zm^pm/Ip-r^tM] 

in the limit of low densities, because equations for the vapor 
pressures, P a (T), and for the saturated vapor densities, 
p a (T), generally are independent. From equations which we 
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used for propane, for example [6], we find that at p = 0, 
(hence T p (p) = 0), the value of Z a {p) may range from to 
00 . Even if Z a (p) — » 1 at low densities, there remains the 
problem of loss of significant digits in the difference [Z a (p) 
~ I]- 

3. Saturated Vapor Densities 

A solution for the problem of eq (2) is to describe saturated 
vapor densities in terms of their compressibility factors 
(utilizing the vapor-pressure equation) in such a way that the 
integrand in eq (2) will extrapolate to a reasonable value as 

The following formulation utilizes temperature, T, as 
independent variable. If, for a given vapor density, however, 
the temperature has been found from this formulation by 
iteration, then T -» T & (p), and the compressibility factor 
becomes a function of density, Z .(7 T ) — > Z o .[7 1 .(p) ], as 
needed for eq (2), above. 

Define the constant A = (Z c — l), where Z c = 
P c /(p c -R* -T c ) is value of the compressibility factor at the 
critical point, and the reduced variables, 



11(f) - P„(T)/P e , x(T) m T/Tc u(T) m l 
when the equation for saturated vapor densities, 

is given by 

Z a (T) = l+A -n-x-* -f(x), 

5 

fix) = 1 +A 1 -u € + ^ J A i >u i -\ 



X, 



(3) 

(4) 
(5) 



At T = 0, f{x) is finite, and Z a (T) = 1, due to the 
exponential dependence of the vapor pressure, P a (T) ~ 
exp(-l/r). 

To give substance to these analytical results, we present 
the vapor-pressure equation for ethylene in appendix A, 
taken from an unpublished report. Figure 1 shows behavior 
of/(V), eq (5), for which the constants are 



e = 0.36 



A 3 = 0.6955 385 



A r = -0.8648 073 A 4 = -3.0648 51 

A 2 = 0.8931 234 A 5 = 5.8633 40 

The rms relative deviation of the ethylene saturated vapor 
densities is 0.07 percent for a total of 29 points, including 
experimental data from 233 through 282 K, and data from 
120 through 220 K derived via the vapor-pressure and virial 
equations [7]. 



f(x) 




.2 .4 .6 

X 

FIGURE 1. Behavior off(x), eq (5), for ethylene. 



4. Conclusions 

By use of eqs (3) and (4) for substitutions in eq (2) at 
densities p < p c , one obtains 



A£ 



= f[(z c 

Jo 



l)-(Z„(p)/Z e ) 



■RT e -f(x) + ...]• dp, (2a) 



in which the integrand is finite in the limit p = 0, {T a (p) = 
0). 

Because the equation of state is nonanalytic, the integra- 
tions of eqs (2) or (2a) must be performed numerically. In 
published works [4,5,6] we have not called upon the equation 
of state at densities so very low as to encounter the difficulty 
described and resolved in this report. At very low densities, 
the saturated vapor temperature diminishes roughly only as 
the logarithm of density, 1/T ~ ln(l/p). 

As it is improbable that we ever shall be privileged to 
compare thermodynamic properties from eq (2) with those 
from eq (2a), our objective in this report has been to describe 
the problem, and to give a solution for consideration in all 
future work. 

5. Appendix A. Ethylene Vapor-Pressure 
Equation 

For use in eq (3), above, our unpublished vapor-pressure 
equation for ethylene, from the triple- to the critical-point, 
for P in bar, is 

ln(P) = a + b-u + c-x + d-x 2 + 

e-x 3 +/•*•(! - x)\ (1A) 



326 



with arguments x = (T - T t )/(T C - T t ), u = (1 - T t /T)/(\ 
- T t /T c ), where e = 1.30, and 



a = -6.717 5716 
b = 13.897 6040 



d = 4.553 6920 
e = -1.132 0564 



c = -6.681 2125 /= 0.623 6504 

For 38 data points from 104.9 through 282.3 K, the rms of 
relative pressure deviations is 0.04 percent [7], 
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